A Nonperturbative Approach to One-Particle Green's Function 
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A nonperturbative method to obtain on- and off-site one-particle Green's function is introduced 
and applied to noninteracting Hubbard model with next nearest neighbor hopping and interacting 
Hubbard model in large dimensions, for example. The former gives some lessons on the method 
and shows the advantage of the method compared with continued fraction formalism. The latter 
is treated by selecting important dynamic processes contributing to the Green's function when 
correlation is strong. We consider the model in the Bethe lattice with large connectivity. The 
dynamic processes describing on-site spin fluctuation clearly shows metal-insulator transition in the 
paramagnetic ground state at half filling. 
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One of recent interesting subjects in condensed mat- 
ter physics is the phenomena of the strongly correlated 
systems. A typical characteristic of the strongly corre- 
lated system is the hardness of getting reasonable solu- 
tion in terms of the traditional perturbative method. The 
Hubbard model, for instance, is one of typical strongly 
correlated systems and is recently revisited with great 
interest as a possible model for high temperature super- 
conductor. 1^,1^ This model has not been solved exactly 
except some of statics in one dimension via Bethe ansatz 
even though it looks simple. Recently, however, the stat- 
ics 1^ and dynamics |^-^ of the Hubbard model have 
been treated exactly in large dimensions, where mean- 
field approximation is valid, with help of numerical tech- 
niques. These exact solutions show interesting behaviors 
such as band collapsing and quasiparticle mass enhance- 
ment near Mott- Hubbard transition, which usual pertur- 
bative method cannot show. In finite dimensions, how- 
ever, only numerical methods such as exact diagonaliza- 
tion and quantum Monte Carlo method have shown some 
meaningful results |^ even though they have severe limi- 
tations. Therefore, an effective nonperturbative method 
is badly required. 

This Letter reports a novel method to calculate one- 
particle Green's function nonpcrturbatively. The Green's 
function is usually obtained by summing Fcynman di- 
agrams ||l0| or truncating equation of motion for the 
Green's function. The former is well-known and solid 
in its foundation, but it is perturbative in principle and 
the latter suggested by Zubarev needs truncation at 
very low level. The famous works of Hubbard jl^ have 
been done using the latter method. But the results are 
only partly acceptable. The method we present here is 
very simple in its form and has never been used before as 
far as we know. We test the effectiveness of the method 
by solving finite-size noninteracting Hubbard model and 
interacting Hubbard model in large dimensions. 

In this work, we treat the retarded Green's func- 
tion which is related to other types of Green's func- 
tions such as time-ordered and thermal Green's function 



and to physical quantities such as single-particle den- 
sity of states (DOS) and internal energy directly. The 
well-known one-particle retarded Green's function in real 
space is 



G%{t) = -temc,.{t),cl}) 



(1) 



and its Fourier transform is 



G%{cj + ^e) = ^-J^ {{c,At), cLl) exp^'^* dt, (2) 

where the superscript R denotes the notation of the re- 
tarded Green's function. Let's notice that the main part 
of Eq. (1) is the projection of the time evolution Cj„{t) 
onto Cka in the Liouville space whose inner product is 
defined as 



(/,g)^({/,gt}) 



(3) 



where / and g are vectors of the Liouville space and the 
angular brackets mean ensemble average, and Eq. (2) 
is the Laplace transform of {{cjait), c\^}) . Thus our in- 
terest to get the Green's function (1) or (2) is in the 
projection {{cjait), clj) . ^ 

Let us assume that we have a complete set of or- 
thogonal bases {e^|^ = 1,2, •••,00} where Ci 



which spans the Liouville space describing the dynam- 
ics of Cjcr(i). Thus the time evolution ei(t) = Cja-{t) can 
be described in terms of these bases. 



ei(i) 



E 



(4) 



Here Ai{t) = (ei(t), ei)/(ei, d) = {{cja{t),cl^}) is just 
the projection we are interested in Eqs. (1) and (2). To 
get Ai{t), we apply the operator identity d/dt — iL to 
Eq. (|^), where L is the Liouville operator and we set 
h = 1. Then one obtains 



^ Ay{t)iLe^ 



(5) 
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Taking inner product in both sides of Eq. (^) with e^, 
one gets a system of hnear equations for A^, (t) such that 



(6) 



Laplace 



for /i — 1, 2, • • • , oo, where M^j^ 

transform of Eq. (^) is written as zAi{z) — 
1 = ^^if^i'(^) for ^ = 1 and zA^{z) = 

^^j^ M^j^Ai^(z) for /i > 2. Use of the boundary con- 
ditions = 0) = 1 and A^{t = 0) = for > 2 
has been made. These equations can be represented in a 
matrix form 



(zI-M)Ac = Bc, 



(7) 



where I is unit matrix and Ac and Be are column matri- 
ces with elements Ai{z),A2{z), ^3(2:), • • • and 1, 0, 0, • • •, 
respectively. Using Cramer's rule, one gets 



A,iz) 



cofactor of (2I — M)ii 
det{zl - M) 

({c,<xW,4})exp-^*dt 



(8) 



and 



A^z) 



cofactor of (zl — M)i 
det{zl - M) 



(9) 



{{e^,el}) Jo 



({c,aW,4}>exp-"*di 



for > 2, where numerators are cof actors of the elements 
in det{zI-M). 

It is well-known that the single-particle DOS P(y{Lo) is 
given by the one-particle retarded Green's function at the 
same site |12], i.e. 



- lim Re^i(z) 



(10) 

(11) 



On the other hand, the off-site retarded Green's function 
Gfki^ + is given by A^{z) of Eq. (9) when 
Therefore, our task to get one-particle retarded Green's 
function is to construct orthogonal bases describing the 
dynamics of the fermion operator, Cja(t], first and calcu- 
late the matrix elements M^^ of Eq. (0) . 

This is a right place to give a remark on this formalism. 
One may consider the operators {iL^ei as a complete set 
of linearly independent bases and construct orthogonal 
bases in terms of Gram-Schmidt process. Then we have 
the following recurrence relation |jl4|--|l6| 



(12) 



where eo = ^^ol^ 



and Ai, 



_ (e^ + i,e^+i) 



when the inner product satisfies the relation (iLe^, c^) = 
~{e^,iLe^) which implies the time-reversal symmetry of 
the inner product. Therefore, when inner product con- 
tains equilibrium average such as Eq. (^ , this relation is 
always satisfied. 

The recurrence relation (|l^) makes the matrix ele- 
ments Mfj^i, vanish for — ^| > 2, since inner products 
{iLei,,e^) vanish, and this reduces the determinants of 
Eq. (8) to tridiagonal forms. Then Ai{z) can be written 
as a well-known infinite continued fraction |p"^-jl6|] 



Aiiz) = 



ai 



(13) 



-012+- 



However, it is practically hard to construct orthogonal 
bases using {iL^ei and impossible to get off-site Green's 
function. More recommendable way is to find linearly 
independent bases intuitively with help of {iL)^ei which 
generates all kinds of linearly independent operators. 

The author and Kee ||l^ have studied the Hubbard 
model using recurrence relation ( p^ and continued frac- 
tion (13) in terms of the bases obtained by large-J7/t* 



expansion. Even though the bases were orthogonal, they 
did not satisfy the recurrence relation for whole region of 
U /i*. Therefore, the continued fraction formalism cannot 
describe the physics of intermediate and low U/t^ regime. 
For this reason, the experiment for insulating V2O3 |lj] 
has been explained very well. |l^ Similar method has 
been apphed to doped Hubbard model with large-J7/i, 
and shown characteristic features of La- and Nd-based 
doped cuprates quite well. (20) To describe the transition 
regime well one must use present formalism instead of 
continued fraction (p^). 

We present here a tutorial example of our formalism us- 
ing noninteracting Hubbard model on 4-site lattice with 
next nearest neighbor hopping, i.e. 



Ho = -t{c\c2 + 4^3 + 4*^4 + c\ci + h.c.) 
+ t'{c\cz + 4c4 + h.c). 



(14) 



We first consider t' = Q case. One naturally selects 
ci,C2,C3,C4 as 4-orthogonal bases 61,62,63,64. Using 
these bases Eq.(8) gives Ai{z) = (z^ + 2t^)/[z{z'^ + U"^)] 
and p(cj) = \5{uj) + \5{uj - 2t) + \5{uj + 2t). On the 
other hand, ^2(2) and ^3(2) give off-site Green's func- 
tions such as 0^2(2) = ^i^ipiF and Gf^{z) = ^(^^It^^) ' 
respectively. One can also easily construct orthogonal 
bases using [iLYci, i.e., Eq. (p^ . If we choose 61 — ci, 
then Eq. ( p^ ) gives 62 — it{c2 + 64), 63 = —2t^C3, and 
64 = 0. (iL^ci for v > 3 are no longer linear indepen- 
dent vectors. Eq.(|l^) and Eq.(|l^) give rise to the same 
Ai{z) and p{uj) as above. However, this does not give 
off-site Green's function. 
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Including next nearest hopping discriminates two 
methods. The first approach raises no difficulty in 
treating next nearest hopping, while the second raises 
big complexity in constructing orthogonal bases. The 
first method easily gives Ai{z) — [z^ + {2t^ + t'^)z + 
2itH']/[z^ + (4^2 + 2t'2)z2 + mH'z - UH'^ + 1'^] from Eq. 
(8). 

This trivial example gives us some lessons about this 
method. First, the bases represent possible independent 
ways of annihilating fermion at a particular site. Second, 
the bases obtained intuitively also give the same result as 
one obtained by using recurrence relation and continued 
fraction, while this does not give off-site Green's function 
and meets great complexity when next nearest neighbor 
hopping is included. Third, one can figure out the mean- 
ing of bases, e.g. t^c^ means annihilation at site 1 after 
twice hopping from site 3. Therefore, one can construct 
bases describing Cj{t) by considering the roles of bases 
intuitively. It is safe to be guided by Eq. ( |l2|) when one 
constructs independent bases intuitively. 

As a meaningful example of constructing bases by 
considering their roles, we treat the half-filled Hubbard 
model in large dimensions. In describing Cja{t), where a 
denotes up-spin, one must consider all possible processes 
while an up-spin electron at site fc, which represents all 
lattice sites including j, hops, arrives at site j , and anni- 
hilates. For the strongly correlated case [U > t), where 
U and t are on-site Coulomb interaction and hopping 
integral, respectively, the down-spin number fluctuation 
(on-site spin fluctuation) when up-spin exists at the same 
site is the most effective part in a basis representing the 
process of the same order of energy, because it yields U 
for each fluctuation. 

We consider the Bethe lattice with arbitrarily large 
connectivity, which is rather easy since the mean-field ap- 
proximation is valid. We show the metal-insulator tran- 
sition in the paramagnetic ground state for the half-filled 
Hubbard model by calculating the single-particle density 
of states ([n] ) . This is the well-known phenomenon stud- 
ied by Hubbard |]l2| and Brinkman and Rice long 
time ago. Recently the infinite dimensional model has 
been studied exactly by several authors. |^-|| 

The interacting Hubbard model is written as 



U 



(15) 



We consider only nearest neighbor hopping. To construct 
bases, we choose ei = Cj„ as before and consider down- 
spin fluctuation at each site where up-spin electron hops. 
We collect the dominant processes describing down-spin 
number fluctuation and neglect hopping dominant pro- 
cesses such as hopping without down-spin number fluc- 
tuation. The bases are written as 



(16) 



X {6n-j,^a + SUq^^^) ■ ■ ■ {5np,^a + 5nq,^a) SUq^^crCq.a, (17) 

for v > 1, where Srij^^^ — nj^^a- — (n^ i> above the 
overbrace denotes the number of the same symbols, and 
the prime means that all site indices under summation 
sign are different each other. The orthogonality of above 
bases is retained by site-difference and {Srii^^a) = 0. The 
property 6nf — {Sui^^^ = 1/4 for half-filhng is used. 

The bases of Eqs. (|l|) and (0) do not satisfy the 
recurrence relation ( |l2[ ) even though they are orthogo- 
nal each other. Therefore, one cannot use the continued 
fraction formalism ([l3| ) . One must use more unrestricted 
formalism (8) to calculate one-particle on-site Green's 
function. The matrix elements Af^t^ are obtained as fol- 
lows in terms of the definition of inner product (|^) and 
the mean-field approximation neglecting spatial correla- 
tions, e.g. {nj^ank,a) = {nj^a){nk,a) for j ^ k, which is 
valid in the limit of large connectivity: 



M2v-l^2v+2' 



iU 



M. 



C/2 



2i/-l,2iy 



^'2' M2u+l,2v — 2' 



-, M2u+2,2v~l — — 



C/2 



(18) 



for V >\, where we scale t as t^j where q is the co- 
ordination number of Bethe lattice, to make the kinetic 
energy finite for large coordination. All other matrix ele- 
ments vanish because electron leaving j site cannot come 
back to the original site via nearest neighbor hopping. 
In calculating the matrix elements, use of the approxi- 
mation g — 1 « q in addition to mean-field approxima- 
tion has been made, which is good enough for a system 
with large connectivity. Once matrix elements Af^^ are 
obtained, the single-particle DOS can be given by Eq. 

(Ill)- 

Fig. 1 shows the results of the single-particle DOS 
for various C//t,. The rank of the matrix M is 101. The 
metal-insulator transition is clearly shown at \Jc = 2t* by 
showing band collapsing from insulating side and quasi- 
particle mass enhancement via disappearance of Fermi 
liquid quasiparticle band from metallic side. 

As a conclusion, we introduced a nonperturbative 
method calculating the one-particle Green's function and 
applied it to a simple tutorial model and the half-filled 
Hubbard model on a Bethe lattice with large connec- 
tivity. The former provides some lessons on this new 
method and shows advantages compared with contin- 
ued fraction formalism, while the latter shows interesting 
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Mott-Hubbard transition by collecting orthogonal bases 
describing important dynamical processes of strong cor- 
relation. 
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The most crucial part of this method is the construc- 
tion of orthogonal bases which describe the various situ- 
ation of electron annihilation processes. These processes 
are obviously very complicate. However, as an approxi- 
mation, one can select some important processes which 
play a major role in a particular parameter regime. The 
second example shows this well. Therefore, this formal- 
ism may provide a possibility to treat the strongly cor- 
related systems in lower dimensions, which is a great 
challenging problem in current condensed matter physics. 
More applications of this formalism will follow. 

This work has been supported by BSRI 97-2420 of the 
Ministry of Education and the KOSEF through the SRC 
program of SNU-CTP. 




1.5 -1 -0.5 0.5 1 1.5 0) 




0.3 



0.2 



0.1 



3 CO 



(d) 






/ \ 
/ \ 

\ 

/ 

\ 


U=3.0 






\ i 

V y 




1 



3 (0 



FIG. 1. The single-particle DOS for various U. We set 

= 1. The (5-function width parameter is chosen e = 0.04 
except in the quasiparticle band in (b) where e oc — ep)]^ ■ 
(b) is magnified to show the band clearly. 
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